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Abstract 



We use the results of direct and indirect searches for weakly interacting mas- 
^ . sive particles (WIMPs) to obtain bounds on various electro-magnetic form factors of 

WIMPs. The limits on the recoil signal in underground dark matter detection experi- 
ments, with standard assumptions about the density and r.m.s. velocity of WIMPs in 
the halo, can be translated into the following model-independent bounds on the mag- 
netic dipole, electric dipole, quadrupole and anapole moments, the charge radius, and 
the polarizability of WIMPs: n D /u n < 1.4 10~ 4 , d < 5 10" 21 ecm, Q < 6 10" 7 efm 2 , 
a/nn < 3 1CT 2 fm, r 2 D < 9 1(T 7 fm 2 , and x < 5 1(T 8 fm 3 for a WIMP mass of lOOGeV. 
The limits on fluxes of highly energetic neutrinos produced in the annihilation of 
WIMPs in the center of the earth and the sun lead to somewhat stronger, but more 
model dependent bounds. 



1 Introduction 



Various observational facts, notably the dynamics of spiral galaxies, clusters and super- 
clusters, and their theoretical explanations in modern cosmology and astrophysics can be 
reconciled only by postulating the existence of the non-luminous matter at amounts ex- 
ceeding the total matter density provided by baryons. Thus, if gravity does not undergo 
some drastic changes at distances larger than few kpc (a highly improbable option from 
various points of view), the explanation of the rotational curves of spiral galaxies requires 
the existence of halos, composed of dark matter objects. 

One of the most intriguing options from the point of view of particle physics is the pos- 
sibility that dark matter is composed of massive particles of yet unknown identity. On the 
one hand these (presumably) weakly interacting massive particles (WIMPs) are the subject 
of intense experimental searches at numerous underground detectors |IJ . On the other hand, 
some of the extensions of the Standard Model naturally incorporate such stable particles. 
For example, in supersymmetric extensions of SM the lightest supersymmetric particle can 
be stable and neutral, thus providing a viable dark matter candidate. Extensive theoretical 
studies of stable SUSY particles, in particular neutralinos, have been performed over the 
last two decades U for wide ranging values of the supersymmetric masses and couplings. 
We believe, however, that unless supersymmetry is experimentally verified and the stabil- 
ity of the lightest superpartner is established in collider experiments, other possibilities for 
WIMPs should not be discarded. 

In this letter we adopt a completely different, yet fully justified, phenomenological ap- 
proach to the physics of WIMPs. We do not make any particular assumptions about the 
origin and possible identity of these particles, apart from the assumption that their masses 
are larger than a few GeV. We then assess what can be inferred about the properties of 
WIMPs directly from experimental limits implied by various searches. An additional incen- 
tive to address WIMP physics in a model independent way stems from recently discovered 
problems with subgalactic structure formation within the non-interacting cold dark mat- 
ter scenario A "generalized" WIMP with sufficient self- interaction may avoid these 
problems Q. 

Among important low-energy properties of WIMPs are their electromagnetic form fac- 
tors. It has been known for a long time that the possibility of charged WIMPs is strongly 
disfavored and stringent limits exist in the case of a fractional charge ||. Here we 
assume that WIMPs are neutral and discuss other electro-magnetic form factors such as 
magnetic and electric dipole moments, anapole moment and polarizability. We put strong 
constraints on the possible electromagnetic form factors of WIMPs using the result of di- 
rect experimental searches aimed at the detection of the WIMP-nucleus recoil signal. These 
constraints depend on the number density and average velocities of WIMPs in the halo. 
The limits on fluxes of energetic neutrinos, originating from the annihilation of WIMPs in 
the center of the Sun and the Earth, can also be used to obtain indirect constraints on 
the form factors. These constraints are less certain, as they require additional assumptions 
about annihilation rates and neutrino branching ratios. 
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This paper is organized as follows. The calculations of the WIMP-nuclei cross sections 
are performed in the next section. In section 3 we present direct limits on the form factors, 
using the results of various underground experiments. In section 4 we obtain indirect limits 
on the form factors, using limits on neutrino fluxes. Our conclusions are reserved for section 
5. 

2 Electromagnetic form factors of WIMPs 

The interaction of any compact object and a slowly varying electro-magnetic field can be 
parametrized via a set of electro-magnetic form factors. In our case, the object is a WIMP 
with spin S, and its interactions, linear in external fields B and E, can be written as 

/r = - /J B.|-dE.|-<d-| (i) 

-45(2^1) 1^ + S jSi - ^SiS + 1)] (Q°-E, + «A Sj ) + .... 

Here /i and d are magnetic and electric dipole moments, and M and Q are magnetic and 
electric quadrupole moments. The anapole moment a is the form factor which describes 
the contact interaction with the external current density j 0. 

If the spin of the WIMP is zero, these moments do not exist and the interaction with 
electromagnetic field is given by the charge radius rjj of the WIMP and the polarizabilities 

H = —erl-^-Ei - ^ X eE 2 - X - Xb B 2 - Xeb E ■ B + ... (2) 

If time reversal is a good symmetry, then d, M and xeb necessarily vanish. Conserved 
parity also forbids the existence of a. Here we do not assume that these symmetries are 
preserved a priori, and we allow all of the form factors to exist. 

In what follows we calculate the cross section of the elastic scattering of a WIMP off 
a nucleus due to the existence of the electromagnetic form factors. We start with the 
scattering due to the magnetic moment which is arguably the most interesting case. 

The calculation can be done in the non-relativistic limit, as v/c ~ 10 -3 on average for 
WIMPs in the halo. We assume that the possible size of the form factors is of the order of, 
or smaller than, the respective characteristic values for nuclei, fi < fi n , Q < e fm 3 , and so 
on. It can then easily be checked that this range of velocities also ensures the applicability 
of the Born approximation. In our units, h — 1, c = 1, e 2 = a — 1/137 and /z n = e/2m p . 

The cross-section for the scattering of a Wimp through the interaction of its magnetic 
moment with the magnetic moment of the nucleus is 

„=f«S WD)AA (£±M±i), (3) 
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where Mr(N,D) is the reduced mass of the WIMP-nucleus system, \in and /xjv are the 
magnetic moments of the WIMP and the nucleus, respectively, and / is the spin of the 
nucleus. The cross-section for the interaction between the electric dipole moment of the 
WIMP and the charge of the nucleus is 

a = 8vrZ 2 ^V^y^ln^. (4) 

\e) \vj 36 q min 

It is logarithmically divergent at large distances. In practice, q max can be chosen to be of the 
order of Mr(N, D)v and q m in is determined either by the inverse of the impact parameters at 
which the nuclear charge is screened by the electrons or by the lowest momentum transfer 
that can be detected in the experiment, whichever is larger. The cross-section for the 
interaction between the nuclear electric current and the anapole of the WIMP is given by 

where \x n is the nuclear magneton, and m p is the proton mass. The contribution of the 
nuclear spin current is negligibly small. The cross-section for scattering through the inter- 
action between the electric quadrupole moment of the WIMP and the electric charge of the 
nucleus takes the form 

.^^JlZVfgf^f^. (6) 
9 KV ' ; Ve/ 105(25-1) v ; 

In Eqs. the normalization of the WIMP spin dependence is such that it is equal to 

1 for the lowest spin for which these moments may exist, (S + 1)/3S = 1 for S = 1/2 and 
(S + 1)(2S + 3)/10S(2S - 1) = 1 for S = 1. The cross-section for scattering of a WIMP 
because of its charge radius in the electric field of the nucleus is given by 

a = ^MUN,D)ZW D . (7) 

Finally, the cross-section for the scattering of a WIMP through its polarizability in the 
electric field of the nucleus, is 

a„ 1 -^MUN,D)Zw4- (8) 

This cross-section diverges for a point-like nucleus. We therefore took the charge distribu- 
tion of the nucleus to be homogeneous in a sphere with radius ro ~ 1.2 fm \^A. 

Scattering due to the magnetic quadrupole moments and magnetic and mixed polariz- 
abilities, Xeb, "was also considered. These effects turned out to be additionally suppressed 
and do not lead to interesting limits. We therefore do not quote them here. 
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3 Direct constraints 



We illustrate the bounds on the electro-magnetic form factors of WIMPs which are implied 
by the limit on the recoil signal with results of the Nal detector of the DAMA experiment 
||, and Ge detector of the CDMS collaboration ||. Other experiments aimed at the 
direct detection of dark matter can be used as well |L0], [□], to get limits which would 
typically be a few times weaker (see also Ref. Q for more complete listing of experiments). 
Traditionally, the results of the direct searches are presented in the form of a normalized 
bound on the "WIMP-proton" cross-section in the spin dependent case, and a normalized 
"WIMP-nucleon" cross-section in the spin independent case. Such normalized bounds are 
useful to compare results for different types of nuclei, to compare results from different 
experiments, and to compare theoretical predictions with experimental data, but for our 
purpose of determining bounds on the WIMP electro-magnetic form factors, it is necessary 
to reverse this normalization. The bounds, quoted in ||, [TIJ are, in some sense, specifically 
formulated for neutralino-nucleon scattering, which necessarily implies a nuclear theory 
input when switching from nucleus-neutralino to nucleon-neutralino description. It is clear 
that in the case of the scattering due to the WIMP electro-magnetic form factors nuclear 
theory is not involved, and much of the nuclear physics uncertainty is not there. Thus, in 
order to obtain bounds on the form factors, we recover the nucleus- WIMP cross section 
from the nucleon-WIMP by the use of the equation || 

M 2 R (N,D) I(N) 
M 2 R (p,D) I(p,n) 

Here p and n stand for proton and neutron, respectively, and the ratio of the reduced 
masses takes into account a kinematic factor. In the spin dependent case, I(N)/I(p) is a 
spin factor, which takes the values I(N)/I(p) = 0.009 for Iodine, and I(N)/I(p) = 0.055 
for Sodium |10|, [13]. In the spin independent case, I(N)/I(n) = (A — Z) 2 for Ref. |10| and 



° = ^2 7Z n\ Tl \ a P,n- I 9 ) 



I(N)/I(n) = A 2 for Ref. |J. For simplicity, we set all nuclear form factors equal to one, 
which is a better approximation for Sodium than for IodineJ] A more accurate analysis 
could be done by the experimental collaborations using their original data. The ensuing 
bounds on the WIMP electro-magnetic form factors are listed in Table (0). 

The bounds presented in this table are obtained with the standard assumptions that 
the mass density of WIMPs in our halo is 0.3 Gev/cm 3 and that their velocities are of 
the order of 200 km/sec. The data obtained with germanium detectors are sensitive to all 
electro-magnetic form factors discussed here except the magnetic dipole moment, as the 
most abundant germanium isotope has no spin. Let us assume that the recently reported 
annual modulation of the DAMA signal, consistent with mo ~ 50 GeV [^4| , and the absence 
of the signal, resulting in strong limits on spin-independent cross section, quoted by CDMS 
||, are both correct. Then the presence of a WIMP magnetic moment at the level of 
few x 10~ 5 - 1 x 10 -4 compared to the size of the nuclear magneton can reconcile the two 
experimental results. 

1 For a more rigorous treatment electric and magnetic form factors of nuclei should be included, for many 
of which the experimental results are available. 
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Na I 


Ge 


y/(S + l)/3S 


1.4 10" 4 




\d/e\J(S+l)/3S 


8 l(T 21 cm 


5 l(T 21 cm 


\Q/e\ yJ(S + 1) (2S + 3) /10S (2S - 1) 


1 10- 6 fm 2 


6 10~ 7 fm 2 


\ a IHn\ \/('S' + i) z^s 


4 10" 2 fm 


3 l(T 2 fm 




1.4 10- 6 fm 2 


9 lO^fm 2 




5 10- 8 fm 3 


7 10- 8 fm 3 



Table 1: Bounds on the electro-magnetic form factors for WIMPs with a mass of 100 GeV, 
based on data from the DAMA collaboration presented in Ref. || and from the CDMS 
collaboration presented in Ref. [[|. 

It is instructive to consider how these bounds are relaxed in the scaling regime, when the 
WIMP mass is larger than nuclear masses and their abundance is kept as a free parameter. 
To this end we introduce the parameter Sd, characterizing the fraction of WIMPs of species 
D in the halo dark matter density (5d = 1 in Table 1). The counting rate in a detector is 
proportional to the number density of the WIMP particles and therefore should scale with 
Md and 5d according to the following formula: 

R n , d ($d, M d ) = R° N 5 D Ij^J , (10) 

where is the lowest mass of WIMPs at which the scaling sets in and -R^v i s ^ ne counting 
rate, corresponding to M D = Mp and 5p = 1. It is easy to see from the figures in Refs.|| [| 
that M D = 100 GeV is close to the mass where the scaling behaviour begins || |], and 
therefore we can simply rescale the bounds given in Table 1 for an arbitrary Sd and arbitrary 
and heavy Md, by multiplying them by ^/jcMo/lOOGeV. 

4 Indirect constraints 

The elastic scattering of WIMPs in the solar (and earth) media leads to gravitational 
trapping. Due to the scattering, WIMPs dissipate energy and eventually they accumulate 



5 



in the center of the sun (and earth). Their density grows until the accumulation rate 
is counter-balanced by the pair annihilation rate, and an equilibrium is reached. In the 
process of annihilation, among other SM products, energetic neutrinos can be produced. 
The interaction of these neutrinos in the earth in turn leads to the production of energetic 
muons, which can be detected with neutrino telescopes. 

Thus limits on the flux of energetic neutrinos can be used to set bounds on the electro- 
magnetic form factors of WIMPs. For WIMPS interacting through their electro-magnetic 
form factors, the trapping rate in the center of the sun and earth is determined from the 
elastic WIMP-nucleus cross sections that we calculated in Eqs.(|3])-(|8]). 



We use the results of Ref. |[L5|| , which follows the original treatment in Refs. and 
estimate the neutrino flux at the surface of the earth due to WIMP annihilation in the sun 

as 

(p uQ ~ (hmcm~ 2 s~ l )N eff 5 D a p ^— —. (11) 

M D 

In this formula N e ff is the average number of neutrinos per annihilation event and <7 Pj 36 
is the WIMP-proton elastic cross section in units of 10~ 36 cm 2 . It is assumed that the 
equilibrium number of WIMPs in the sun has been reached, so that the annihilation and 
capture rates are equal. 

The experimental limit on the flux of energetic upward muons (for example obtained by 
the Kamioka, MACRO and Baksan experiments, [|T7|, [R| [T§), is fa < 1.4 x 10" 14 cui'V 1 , 
and the probability that a neutrino directed towards the detector produces a muon at the 
detector is P(lOOGeV) ~ 1CT 7 p0|| . Using Eq . (|TT|) , we deduce the following indirect limits 



on the electro- magnetic form factors of WIMPs with mp = 100 GeV, 



\fi D /^n\\/(S + l)/3S < 5HT 6 
\d/e\ J{S + 1)/3S < 2 10~ 22 cm 



\Q/e\yJ(S + l){2S + 3)/10S(2S-l) < HT 8 fm 2 (12) 

\a/t*n\ y/(S + l)/3S < 10~ 2 fm 

r 2 D < 1.4 10~ 8 fm 2 
|x| < 5.6 10~ 7 fm 3 , 

where N e ff is taken to be 0.1 and djj = 1. Note that for higher WIMP masses the absorption 
of neutrinos in the core of the sun becomes important ||21|| . 

If the WIMP-nucleus elastic scattering cross-sections in the earth are comparable to 
the WIMP-proton cross-section, than the neutrino signal from the center of the earth is 



approximately four orders of magnitude weaker than the signal from the sun |15[ . However, 
the limit on xe from the neutrino flux from the center of the earth is almost one order of 
magnitude better than the limit derived from the neutrino flux from the sun, as the cor- 
responding WIMP-nucleus elastic scattering cross-sections are enhanced by a large factor, 
Z A ~ 5 10 5 , for heavy (Fe) nuclei (see Eq.(|J)) comprising the core of the earth. 
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Although indirect bounds are definitely stronger, they do require a number of additional 
assumptions. In particular, these bounds require that annihilation actually occurs, which in 
practice may not be true. For example, if a WIMP caries a conserved quantum number, and 
there are no "anti-WIMPs" - such as heavy Dirac neutrinos in the absence of antineutrinos 
- the annihilation is forbidden. 



5 Discussion 

We have shown that limits on recoil signals from underground detectors can be translated 
into bounds on the electro-magnetic form factors of WIMPs. Similar bounds can be inferred 
from the limits on energetic neutrino fluxes from the center of the earth and the sun, 
generated by annihilating WIMPs. Our analysis shows that all electromagnetic form factors 
have to be very small when expressed in characteristic nuclear units. In particular, the 
magnetic moment of a WIMP has to be five orders of magnitude smaller than the magnetic 
moment of a proton. 

Unfortunately, all quoted bounds are trivially satisfied for the case of neutralino dark 
matter. Indeed, because it is a spin- 1/2 particle, a neutralino cannot possess a quadrupole 
moment. In addition, because it is a Majorana fermion, a neutralino cannot have a mag- 
netic or electric dipole moment either {x a nvX = 0)- Finally, an anapole moment and a 
polarizability can be generated only radiatively. This results only in a small correction to 
the dominant neutralino-nucleus elastic cross section induced by squark exchange ||. 

A magnetic moment can arise naturally if WIMPs are Dirac fermions. If we assume 
for a moment that the WIMP is a heavy Dirac neutrino, a large magnetic moment can 
be induced due to virtual M^-charged lepton exchange. The value of the induced magnetic 
moment grows quadratically with the size of the Yukawa coupling until the point where 
the Yukawa interaction becomes essentially strong. In latter case a natural estimate for 
the size of the magnetic moment is e/m„, with m„ presumably of the order of 1 TeV or 
heavier. In this case, the interaction through the magnetic moment can be very important 
and, for light nuclei, can even dominate Z-boson exchange. It is interesting to note that 
the electric dipole moment of a WIMP is also very constrained. The limit of 5 10~ 21 e cm 
is stronger than any other EDM constraint, apart from EDMs of electrons, neutrons and 
heavy atoms. The method we used to obtain this limit is in the spirit of the original paper 
by Purcell and Ramsey |22] , where the first bound on the neutron EDM was inferred from 



neutron-nucleus scattering. 

A model independent approach to the physics of WIMPs is partly motivated by recently 
discovered problems within the cold dark matter model. While large-scale structure can be 
described reasonably well, the numerical simulations of galactic substructures appear to be 
in conflict with observational data ||. A recently proposed cure, invoking self-interacting 
dark matter M , is incompatible with the very restrictive neutralino model. This fact gives 
further impetus to our approach to study WIMPs in a manner as model independent as 
possible. The limits on the electro-magnetic form factors obtained in this work show that 
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the amount of self-interaction, which could be induced by these form factors, is not sufficient 
to generate a cross section of WIMP-WIMP scattering at the level of 10 -25 cm 2 as required 
by the self-interacting WIMP scenario. 

The authors would like to thank Keith Olive, Arkady Vainshtein and Viktor Zacek for 
usefull discussions. This work was supported in part by the Department of Energy under 
Grant No. DE-FG-02-94-ER-40823. 
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